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We address generation, propagation and application of multipartite continuous variable entan- 
glement in a noisy environment. In particular, we focus our attention on the multimode entangled 
states achievable by second order nonlinear crystals, i.e. coherent states of SU(m, 1) group, which 
provide a generalization of the twin-beam state of a bipartite system. The full inseparability in the 
ideal case is shown, whereas thresholds for separability are given for the tripartite case in the pres- 
ence of noise. We find that entanglement of tripartite states is robust against thermal noise, both 
in the generation process and during propagation. We then consider coherent states of SU(m, 1) 
as a resource for multipartite distribution of quantum information, and analyze a specific protocol 
for telecloning, proving its optimality in the case of symmetric cloning of pure Gaussian states. We 
show that the proposed protocol also provides the first example of a completely asymmetric 1-nn 
telecloning, and derive explicitly the optimal relation among the different fidelities of the m clones. 
The effect of noise in the various stages of the protocol is taken into account, and the fidelities of the 
clones are analytically obtained as a function of the noise parameters. In turn, this permits the opti- 
mization of the telecloning protocol, including its adaptive modifications to the noisy environment. 
In the optimized scheme the clones' fidelity remains maximal even in the presence of losses (in the 
absence of thermal noise), for propagation times that diverge as the number of modes increases. In 
the optimization procedure the prominent rule played by the location of the entanglement source 
is analyzed in details. Our results indicate that, when only losses are present, telecloning is a more 
effective way to distribute quantum information then direct transmission followed by local cloning. 

PACS numbers: 03.67.Mn, 03.67.Hk 

I. INTRODUCTION 

Entanglement plays a fundamental role in quantum information, being recognized as the essential resource for 
quantum computing, teleportation, and cryptographic protocols. In the framework of quantum information with 
continuous variables (CV) P, Q] the possibility of generating and manipulating entanglement allowed the realization 
of a variety of quantum protocols, such as teleportation, cryptography, dense coding and entanglement swapping. In 
these protocols the source of entanglement is the bipartite twin-beam state of two modes of radiation, usually generated 
by parametric down-conversion in crystals. However, recent experimental progresses Q show that the coherent 
manipulation of entanglement between more then two modes may be achieved with current technology. This opens 
the opportunity to realize a true quantum information network, in which the information can be stored, manipulated 
and distributed among many parties, in a fashion resembling the current classical telecommunication networks. In a 
realistic implementation, entanglement needs to be transmitted along physical channels, such as optical fibers or the 
atmosphere. As a matter of fact, the propagation and the influence of the environment unavoidably lead to degradation 
of entanglement, owing to decoherence effects induced by losses and thermal noise. In this scenario, it is worth to 
study the entanglement properties and the possible applications of multipartite systems in noisy environments, which 
will be the subject of this paper. 

A prominent class of CV states is constituted by the Gaussian states. They can be theoretically characterized in 
a convenient way, and generated and manipulated experimentally in a variety of physical systems. In a quantum 
information setting, entangled Gaussian states provide the basis for the quantum information protocols mentioned 
above. The basic reason for this is that the QED vacuum and radiation states at thermal equilibrium are themselves 
Gaussian states. This observation, in combination with the fact that the evolutions achievable with current technology 
are described by Hamiltonian operators at most bilinear in the fields, accounts for the fact that the states commonly 
produced in labs are Gaussian. Indeed, bilinear evolutions preserve the Gaussian character. As we already mentioned, 
the outmost used source of CV entanglement are the twin-beams, which belong to the class of bipartite Gaussian 
states. In a group-algebraic language, they are the coherent states of the group SU(1, 1), i.e., the states evolved 
from vacuum via a unitary realization of the group. Within the class of Gaussian states, the simplest generalization 
of twin-beams to more then two modes are the coherent states of the group SU(m, 1). Indeed, this states can be 
generated by multimode parametric processes in second order nonlinear crystals, with Hamiltonians that are at most 
bilinear in the fields. In particular, these processes involve m + 1 modes of the field ao,ai,... ,a m , with mode ao 
that interacts through a parametric-amplifier-like Hamiltonian with the other modes, whereas the latter interact one 
each other only via a beam-splitter-like Hamiltonian PI [5] . In the framework of CV quantum information, the first 
proposal to produce such states has been given in Ref. Q , where a half of a two- mode squeezed vacuum state interacts 
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with m vacua via a proper sequence of beam splitters. Other unitary realizations of the algebra of SU(m, 1) have been 
proposed, in optical settings 0, @ as well as with cold atoms @ or optomechanical systems [T^j • ^ n these schemes 
the Hamiltonian of the system, rather then involving a sequence of two-mode interaction, is realized via simultaneous 
multimodc interactions. Experimental realizations of tripartite entanglement in the optical domain have been recently 
reported Q. 

In this work we do not focus on any specific implementation of the SU(m, 1) evolution. Rather, we will analyze 
the entanglement properties of SU(m, 1) coherent states in a unified fashion valid for a generic Hamiltonian of this 
kind. As we will see in Sec. [HI this is allowed by the observation that the coherent states of SU(m, 1) have a common 
structure, which can be conveniently written in the Fock representation of the field In particular, the degradation 
effects of both the thermal background in the generation process and of losses and thermal photons in the propagation 
will be outlined. The robustness of these states against noise will be analyzed in Sec. lIIII where it will be also compared 
with the bipartite case. 

As already mentioned, one of the main results in CV quantum communication concerned the realization of the 
teleportation protocol (for a recent experiment see LO|)- The natural generalization of standard teleportation to 
many parties corresponds to a telecloning protocol Teleportation is based on the coherent states of SU(1,1), 

which provide the shared entangled states supporting the protocol. Thus, in order to implement a multipartite 
version of this protocol, one is naturally led to consider a shared entangled state produced by a generic SU(m, 1) 
interaction. The telecloning protocol will be analyzed in detail in Sec. IIVI As concern cloning with CV, there are 
general results to assess the optimality of n — ► m symmetric cloning of coherent states 13] . Optimal local unitary 
realization of such schemes have been proposed in |l4t Il5| , and an experimental realization of 1 — ► 2 cloning has been 
recently reported fl6| . Concerning telecloning, existing proposals are about optimal 1 — > m symmetric cloning of pure 
Gaussian states, using a particular coherent state of SUf w, 1*) as support [fj. Recently, a proposal which make use 
of partially disembodied transport has also been reported 17] . In view of the realization of a quantum information 
network, one is naturally led to consider the possibility to retrieve different amount of information from different 
clones. This means that one may consider the possibility to produce clones different one from each other, in what 
is called asymmetric cloning. Examples of optimal 1 — * 2 asymmetric cloning are given in Ref. |l5l Il8| . where local 
and non-local realizations are considered. In this work, we will see how the telecloning protocol involving a generic 
coherent state of SU(m, 1) provides the first example of a completely asymmetric 1 — > m cloning of pure Gaussian 
states. In this sense, we provide a generalization of the proposal in Ref. |6( to the asymmetric case. Moreover, we 
found an expression for the maximum fidelity achievable by one clone when the fidelities of the others are fixed to 
prescribed values, thus giving explicitly the trade-off between the qualities of the different clones. 

In Sec. we will analyze the effect of noise in each step of the telecloning protocol. As expected, the presence of 
both thermal noise and losses unavoidably leads to a degradation of the cloning performances. Nevertheless, we will 
show that the protocol can be optimized in order to reduce these degradation effects. In particular, one may optimize 
not only the energy of the entangled support, but also the location of the source of entanglement itself. Remarkably, 
when only losses are considered, this optimization completely cancels the degrading effects of noise on the fidelity of 
the clones. This happens for finite propagation times which, however, diverge as the number of modes increases. 

We conclude the paper with Sec. IVII where the main results will be summarized. 



II. MULTIMODE PARAMETRIC INTERACTIONS: SU(m, 1) COHERENT STATES 



Let us consider the set of bilinear Hamiltonians expressed by 

rn m 

H m = a k a\ + Y l [ k ] a fc°o + h.c. , (1) 



Z<fc=l fe=l 



where [afc,aj] = 0, [a*,, a|] = Sk.i (k,l = 0, ...,m) are independent bosonic modes. A conserved quantity is the 
difference D between the total mean photon number of the mode ao and the remaining modes, in formula 

rn 

D = ^ ajak ~ %ao ■ (2) 

k=l 

The transformations induced by Hamiltonians Q correspond to the unitary representation of the SU(m, 1) algebra 
Therefore, the set of states obtained from the vacuum coincides with the set of SU(m, 1) coherent states i.e. 

|* m ) = exp{-iff m t}|0) =exp j |^ fc 44 -/i.e. j |0> , (3) 
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where (3k are complex numbers, parameterizing the state, which are related to the coupling constants 7^,}' and 7^ in 
Eq. QJ. Upon defining 



C fc = A 



tanh (J27=i \@r 

EZM 2 



\^f m ) in Eq. © can be explicitly written as 



2 ■■■ L 'n 



y/Jjll + n 2 + ... + n m )\ 



\fni\n 2 \...n m \ 



^2n k ;{n}) , 



(4) 



fc=i 



where {n} = {711,713, ••• } 7im}- The sums over n are extended over natural numbers and Z m = 1 — 53fc=i l^fc| 2 is a 
normalization factor. We see that for m = 1 one recover the twin-beam state. Notice that, being interested in the 
entanglement properties and applications of states |^ m ), we can take the Cfc's coefficients as real numbers. In fact 
one can put to zero the possible phases associated to each C% by performing a proper local unitary operation on mode 
ttfe, which in turn does not affect the entanglement of the state. Calculating the expectation values of the number 
operators Nh = (a' k a,k) on the multipartite state |* m ) one may re-express the coefficient in Eq. J3J as follows: 



C k = 



1/2 



-2m — 



1 



(k = 1, 



(5) 



In order to obtain Eq. (JSJ) we have considered Eq. J5J with D — (vacuum input), from which follows that 



fc=l 



and repeatedly used the following identity: 



y x n(n + aV. =al{1 _ x) -i-a 

n=0 



(6) 



(7) 



The case D / will be considered in the next section, in which the effects of thermal background will be taken into 
account. The basic property of states in Eq. Q is their full inseparability, i.e., they are inseparable for any grouping 
of the modes. To prove this statement first notice that, being evolved with a bilinear Hamiltonian from the vacuum, 
the states \^ m ) are pure Gaussian states. They are completely characterized by their covariance matrix <x, whose 
entries are defined as 



(8) 



where {A, B} = AB + BA denotes the anticommutator, R = (qo,Po, ■ ■ ■ , <?m>Pm) T and the position and momentum 
operator are defined as qu = (a& + a k )/\/2 and pk = (<ifc — a k )/y/2. The covariance matrix for the states \^ m ) reads 
as follows: 



/ Af Ai A 2 

A x Mi Si, 2 
A 2 S12 Af 2 



where the entries are given by the following 2x2 matrices (k = 



^m — l,m 



(9) 

m, j = 2, . . . , m and < ? < j) 

(10) 



M"k = (N k + -) 1 A h = y/N h {N + l)P Bij = y/N~N~ 1 , 

with 1 = Diag(l, 1) and P = Diag(l,— 1). Since |^ m ) are pure states, full inseparability can be demonstrated by 
showing that the Wigner function does not factorize for any grouping of the modes, which in turn is ensured by the 
explicit expression of the covariance matrix <j m given above (as soon as Nh ^ 0). 
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III. EFFECT OF NOISE ON THE GENERATION AND PROPAGATION OF SU(m, 1) COHERENT 

STATES 

In view of possible applications of the coherent states of SU(m, 1) to a real quantum communication scenario, it is 
worth to analyze the degrading effects on their entanglement that may arise when generation and propagation in a 
noisy environment is taken into account. Unfortunately, a manageable necessary and sufficient entanglement criterion 
for the general case of a Gaussian multipartite state is still lacking. Thus, in order to study quantitatively the effects 
of noise we must limit ourselves to the case when only three modes are involved (insights for the general m-modc 
case will be given in the following Sections). In fact, u p to three modes the partial transpose criterion introduced in 
jlll I20L l2l| is necessary and sufficient for separability [22j . It says that a Gaussian state described by a covariance 
matrix er is fully inseparable if and only if the matrices uj k — cr — ^Jk {k = 0, 1, 2) are non-positive definite, where 
J fe = A fe JA fc with A = Diag(l, -1, 1, 1, 1, 1), Aj = Diag(l, 1, 1, -1, 1, 1), A 2 = Diag(l, 1, 1, 1, 1,-1) and 



-1 3 
1 3 



(11) 



and 1„ is the nx n identity matrix. This criterion has been applied in Refs. |22ll23| in order to assess the separability 
of the CV tripartite state proposed in Ref. (2^| when thermal noise is taken into account. In Ref. ^(| the entanglement 
properties of a state generated via a SU(2, 1) evolution when one of the modes starts from thermal background has 
been also numerically addressed. 

Let us now analyze if the generation process of states |*P m ) is robust against thermal noise. This means that 
we have to study the separability properties of a state generated by a SU(m, 1) interaction starting from a thermal 
background rather then from the vacuum, in formulae g = e~ %Hmt q v e lffm *, where g and g„ are the density matrix of 
the evolved state and of a thermal state, respectively. We may call these states thermal coherent states of SU(m, 1). 
First notice that, being the thermal state Gaussian, the thermal coherent states will be Gaussian too, and their 
covariance matrix <x mj th may be immediately identified from Eq. @- In fact, in the phase space identified by the 
vector R, every SU(m, 1) evolution will act as a symplectic operation <S on the covariance matrix of the input state, 
i.e., <r out = <S T 0"i n <S. Recalling that the covariance matrix of a thermal state can be written as er th = (2 v + l)cr Vl 
being er v = 1/2 the covariance matrix of vacuum and v the mean thermal photon number, we obtain 

0"m,th = (2v + l)cr m (12) 

Let us now apply the separability criterion recalled above to cr^^h- Concerning the first mode, from an explicit 
calculation of the minimum eigenvalue of matrix luq it follows that 



A™ in = v + (1 + 2u) [N - ^/No(N + 1)J . (13) 
As a consequence, mode oq is separable from the others when 



v > N + ^N Q {N Q + 1) . (14) 
Calculating the characteristic polynomial of matrix lo\ one deals with the following pair of cubic polynomials 

q x {\, N q ,N 1 ,N 2 ,p) = A 3 - 2 [2(1 + N Q ) + ^(3 + 4^)] \ 2 

+ 4 [1 + Nx + 2N 2 + v{A + 4JVi + 6N 2 + v(3 + 4JV ))] A 

- 8^ [1 + JVi + i/(2 + v + 2Ni)] , (15) 

q 2 (A, JV , N u v) = A 3 - 2 [1 + 2N Q + u(3 + AN )] A 2 

+ 4 [Ni + 2i/(l + 7V ) + v 2 {i + 4JV )] A 

-8(l + i/)(i/ 2 -2iVi-2^i) . (16) 

While the first polynomials admits only positive roots, the second one shows a negative root under a certain threshold. 
It is possible to summarize the three separability thresholds of the three modes involved in the following inequalities 



v > N k + ^/N k (N k + 1) . (17) 

If Inequality l|17|) is satisfied for a given k, then mode a k is separable. Clearly, it follows that the state 1*2) evolved 
from vacuum (i.e., v = 0) is fully inseparable, as expected from Section^ Remarkably, Inequality i|17|l is the same 
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as for the twin beam evolved from noise [25|], which means that the entanglement of the thermal coherent states of 
SU(2, 1) is as robust against noise as it is for the case of the thermal coherent states of SU(1, 1). 

Let us now consider the evolution of the state \^ 2 ) in three independent noisy channels characterized by loss rate 
r and thermal photons [i, equal for the three channels. The covariance matrix cr 2 (t) at time t is given by a convex 
combination of the ideal <r 2 (0) [i.e., <t 2 in Eq. JUJ] and of the stationary covariance matrix er oo i2 = (a 1 + 5)^6 

<r 2 (t) = e- r * 0-2 + (1 - e" rt ) 0-00,2 . (18) 

Consider for the moment a pure dissipative environment, namely fi = 0. Applying the separability criterion above to 
er 2 (t), one can show that it describes a fully inseparable state for every time t. In fact, we have that the minimum 
eigenvalue of lvq is given by 



A™ in = 2e- rt 



N Q - y/N (N + 1)1 ■ (19) 



Clearly, A™ 111 is negative at every time t, implying that mode clq is always inseparable from the others. Concerning 
mode Oi, the characteristic polynomial of u>i(t) factorizes into two cubic polynomials: 

q 1 {X,T,N ,N 1 ,N 2 ) = -A 3 + 4 [l + e^'TVo] A 2 

+ 4 [-1 - e - rt (2iVi + 3A 2 - e~ rt iVo)] A + 8e~ r *A^ 2 (l - e~ rt ) , (20a) 
g 2 (A,r,iVo,iVi,iV 2 ) = -A 3 + 2[l + 2e- r *A^ ] A 2 

+ 4 [- e - r *(2Aq + N 2 ) + e" 2r *iV ] A - 8e- 2r *Aq . (20b) 

While the first polynomial has only positive roots, the second one admits a negative root at every time. Due to 
the symmetry of state \^ 2 ) the same observation applies to mode et 2 , hence full inseparability follows. This result 
resembles again the case of the twin beam state in a two-mode channel |20l 1261 ] . In other words, the behavior of the 
coherent states of SU(2, 1) in a pure lossy environment is the same as the behavior of the coherent states of SU(1, 1), 
concerning their entanglement properties. 

When thermal noise is taken into account (fj, ^ 0) separability thresholds arise, which again resembles the two-mode 
channel case. Concerning mode ao, the minimum eigenvalue of matrix u>o(t) is negative when 



/< f lnll + Xi 2 — 2 - J . (21, 



1 



where N to t = Nq + N\ + N 2 . Remarkably, this threshold is exactly the same as the two- mode case [2(j, if one 
consider both of them as a function of the total mean photon number of the TWB and of state |\I/ 2 ) respectively. 
This consideration confirms the robustness of the entanglement of the tripartite state 1*2). Concerning mode ax, the 
characteristic polynomial of cji(t) factorizes again into two cubic polynomials. As above, one of the two have always 
positive roots, while the other one admits a negative root for time t below a certain threshold, in formulae: 

- 8e~ 2r *Ari + 8(e~ r * - l)e~ r *( e - rt 7V - 27Vi - N 2 )fi 

+ 8{e- rt ~ 1) 2 (1 + 2e- rt iVo) j u 2 - 8( e - r * - 1)V < . (22) 

Mode a 2 is subjected to an identical separability threshold, upon the replacement N% ^ N 2 . In Fig. \T\we compare 
the separability thresholds given by Eq. I|21|l and Eq. (|22|) . As it is apparent from the plot, modes a\ and a 2 become 
separable faster then modes ao, hence the threshold for full inseparability of I* 2) is given by Eq. 




FIG. 1: Separability thresholds for modes ao (continuous line) and ai (dashed line) according to Eq. (I21H and Eq. 122H for the 
case of iVi = N2 = N = 1. The behavior of these curves is similar if different values of iV are considered. 
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classical channel 



FIG. 2: Schematic diagram of the telecloning scheme. After the preparation of the state |3?m), a conditional measurement is 
made on the mode do, which corresponds to the joint measurement of the sum- and difference-quadratures on two modes: mode 
ao itself and another reference mode b, which is excited in a pure Gaussian state |£, a), to be teleported and cloned. The result 
z of the measurement is classically sent to the parties who want to prepare approximate clones, where suitable displacement 
operations (see text) on modes a\, . . . , a m are performed. We indicated with v and [i the mean thermal photons in generation 
and propagation, whereas A takes into account for the non unit efficiency in the detection stage. The effective propagation 
times to and r c (see Sec.|VJl are related to the losses during propagation. 

We conclude that the entanglement properties of the coherent states of SU(2, 1) in a noisy environment resembles the 
twin-beam case both in generation and during propagation. This may be relevant for applications, as the robustness 
of twin beam is at the basis of current applications in bipartite CV quantum information. 

IV. TELECLONING 

We now show how the multipartite states |\P m ) introduced in Sec. [H]can be used in a quantum communication 
scenario. In particular we show that states \^ m ) permit to achieve optimal symmetric and asymmetric 1 — ■> m 
telecloning of pure Gaussian states. Optimal symmetric telecloning has been in fact already proposed in Q using a 
shared state produced by a particular SU(m, 1) interaction. Also a protocol performing optimal 1 — > 2 asymmetric 
telecloning of coherent states has been already suggested in Ref. 0] , where the shared state is produced by a suitable 
bilinear Hamiltonian which generates a SU(2, 1) evolution operator. Here we consider the general 1 — ► m telecloning 
of Gaussian pure states in which the shared entanglement is realized by a generic coherent state of SU(m, 1). First 
recall that a single-mode pure Gaussian state can be always written as 

\£,a)=S b (OD b (a)\0) , (23) 

where 5f,(£) = exp {i£(o t ) 2 — |£*o 2 } and A, (a) = exp {ao 1 ' — a*b] are the squeezing and the displacement operator 
respectively, whereas b is the mode to be cloned. We emphasize that our goal is to create m clones of state g m = 
|£, a)(£, a\ in a non-universal fashion, i.e. the information that we clone is encoded only in the coherent amplitude a. 
In other words, we consider the knowledge of the squeezing parameter £ as a part of the protocol, as in the case of 
local cloning of Gaussian pure states j^- The telecloning protocol is schematically depicted in Fig. [3 As a shared 
entangled state we consider the following |2S[ : 

|*m> = S«„(0 ® S ax (0 ® . .. ® S am (£)|* m ) • (24) 

After the preparation of the state |4? m ), a joint measurement is made on modes ao and 6, which corresponds to 
measure the complex photocurrent Z = b + (double- homodyne detection), as in the teleportation protocol. The 
measurement is described by the POVM {H(z)} z£ c, acting on the mode ao, whose elements are given by 

U(z)=TT- 1 D ao (z) <&Dl(z) 

= 7 r- 1 5 Qo (rpa o (^ao(«*)|0)(0|<(a*)<(^ o (r), (25) 
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where z is the measurement outcome, z' = zcoshr + e l8 z* svahr, £ = re 10 and T denotes transposition. The 
probability distribution of the outcomes is given by 



P(z) = Tr . 



1 



|* m )(* m |n(2)(g)i h 

h=l 

\z' + a*\ 2 



exp 



tt(1 + N ) * \ l + N 

being I/! the identity operator acting on mode et^. The conditional state of the remaining modes then reads 

1 



(26) 



0-: 



P{z) 



Trn 



\& m ){* m \n(z)(g)i h 

h = l J 



tS ah (0\C h (z'*+a)){C h (z'*+a)\Sl h (0 



(27) 



h=l 



where \Ch(z'* + a)} denotes a coherent state (of the usual Heisenberg Weyl group) with amplitude Ch(z'* + a). After 
the measurement, the conditional state should be transformed by a further unitary operation, depending on the 
outcome of the measurement. In our case, this is a rn-mode product displacement U z — ®hLi ^^(z). This is a local 
transformation, which generalizes to m modes the procedure already used in the original CV teleportation protocol. 
The overall state is obtained by averaging over the possible outcomes 

Ql...m = / d 2 zP(z) t z . 
JC 

where t z — U z g z U\. Thus, the partial traces Qh = Tri ! ... t h-i,h+i,...,m[Qi...m] read as follows 



Qh = S h (£) 



d 2 z P(z) \a C h + z'* (C fc - 1)) (a C h + z'* (C h - 1) 



slit) 



Upon a changing in the integration variable we obtain the following expression for the clones: 



Qh = s h {Ci 



d 2 w 



1 



irn h 



■ exp 



\w — a\~ 
n h 



\w)(w\ 



where we defined 



n h = ^N Q + 1 - y/N^j 



(28) 



(29) 



(30) 



From expression l|29(l one immediately recognize that the clones are given by thermal states Qth{nh), with mean photon 
number rih, displaced and squeezed by the amounts a and £ respectively, i.e.: 



Qh = S h (0 D h (a) Qth (n h ) D^a) (£) 



(31) 



As a consequence, we see that the protocol acts like a proper covariant Gaussian cloning machine |27j, and that the 
noise introduced by the cloning process is entirely quantified by the thermal photons rih, which in turn depend only 
on the value of the mean photon numbers Nh of the shared state. The fidelity Fh between the h-th clone and the 
initial state does not depend on the latter and is given by 



1 



1 + n h 



(32) 



The expression of the clones in Eq. I|31|l says that they can either be equal or different one to each other, depending 
on the values of the nh's. In other words, a remarkable feature of this scheme is that it is suitable to realize both 
symmetric, when ri2 — ■ ■ ■ = n m — n, and asymmetric cloning, ri2 ^ • ■ • 7^ n rn . This arises as a consequence of 
the possible asymmetry of the state that supports the telecloning. To our knowledge, this is the first example of a 
completely asymmetric 1 — > m cloning machine for continuous variable systems. 
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Concerning the symmetric cloning one has that this scheme saturates the bound given in Ref. 01 > hence ensuring 
the optimality of the protocol. In fact, the minimum added noise for a symmetric 1 — > to doner of coherent states is 
given by n = ^^E^, which in our case can be attained by setting Ni = • • • = N m = N opt , where 

N opt = _ _i (33) 

m(rri — 1) 

It follows that the fidelity is optimal, namely F = to/ (2m — 1). It is not surprising that this result is the same as the 
one obtained in Ref. [|J. In fact, as already mentioned, the latter uses as support a specific SU(m, 1) coherent state, 
generated with a particular interaction built from single mode squeezers and beam-splitters. Our calculation extends 
this result to any SU(m, 1) coherent state used to support the telecloning protocol. 



A. Asymmetric cloning 



Consider now the case of asymmetric cloning. In this case one deals with a true quantum information distributor, in 
which the information encoded in an original state may be distributed asymmetrically between many parties according 
to the particular task one desires to attain. In this scenario, a particularly relevant question concerns the maximum 
fidelity achievable by one party, say Fi, once the fidelities Fj (j — 2, . . . , to) of the other ones are fixed. Thanks to 
Eq. ll-il'l) we see that this is equivalent to the issue of finding the minimum noise n\ introduced by the cloning process 
for fixed rij's (rij ^ 1). The optimization has to be performed under the constrain given by Eq. (|2J, which allows to 
write ni as a function of the rij's and of the total mean photon number A^o (the sums run for j = 2, . . . , to): 



"l 



VNo + 1 - n - ]T (y^o + i - v% 



The minimum noise n 1 ™ 111 is then found setting Nq such that 

l)(m - 1)(to - 2) - 2 y/N + l(m - 1)£,V»£ + + fe^)' 



(No 



1 = 0. 



(34) 



(35) 



ropt _ 



For m = 2 one obtains that the optimal choice for Nq is given by N 
noise n? 11 " allowed by our telecloning protocol for fixed n 2 is given by nj 
Ref. 18] for the fidelities: 



n 2 + l/4n 2 . It follows that the minimum 
mm _ i/4 :Tl2 _ Hence we recover the result of 



pm 



4(1 - F 2 ) 
4-3F 2 



(36) 



Notice that if one requires F 2 = 1 then i*\ = 0, that is no information is left to prepare a non-trivial clone on mode a\. 
We remark that the result in Eq. I|3tj|) shows that the protocol introduced above, besides reaching the optimal bound 
in the symmetric case, is optimal also in the case of asymmetric 1 — > 2 cloning [15|- Coming back to the general case 
we see from Eas. (l34|) and (|35|l that for to > 3 the minimum noise n± is given by 



1 



E 



1 (to-2) 2 
and it is attained for the following optimal choice of Nq 
1 



y/{m-l) [(Z jy /nj) 2 ~ (to - 2) l Z ] n J - (to - 2)] 



N o V t 



(to - 1) 2 (to - 2) 5 



^ (m - 1) [(E 3 >J) 2 - (to - 2)E,^ - (m - 2)] | 



(37) 



(38) 



Substituting Eq. (|37|l in Eq. 132(1 one then obtains the maximum fidelity F™ ax achievable for Fj fixed. Summarizing, 
if one fixes the fidelities Fj (for j = 2, . . . ,to) then the thermal photons nj are given by Eq. I|32|) . which in turn 
individuate the mean photon numbers Nj and Ni of the state that supports the telecloning via Eqs. I|30l) . 13 7|) and 
()38|l . This choice guarantees that the fidelity F± is the maximum achievable with the telecloning protocol described 
above, thus providing the optimal trade-off between the qualities of the different clones. 

As an example, consider the fully asymmetric 1 — > 3 telecloning and fix the couple of fidelities F 2 and F3. Specializing 
the formulae above we have that, choosing the state ^3) such that: 




N 



2(3) 




(39) 
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then the fidelity of the first clone is the maximal allowed by our scheme, in formulae: 



\ 




(40) 



Notice that F^ ax in Eq. igQJ is valid iff the fixed fidelities F 2 and F 3 satisfy the relation F 2 < 4(1 - F 3 )/(4 - 3F 3 ), 
which coincides with the optimal relation given by Eq. (|36|) . In other words, the optimal bound imposed by quantum 
mechanics to 1 — > 1 telecloning is automatically incorporated into the bound (|40|) for 1 — > 3 telecloning of our scheme. 
When F 2 = F 3 = 3/5 (that is, the bound for an optimal symmetric 1^3 doner) we have that F 1 max = 3/5, as one 
may expect from the discussion above concerning the symmetric cloning case. Remarkably, when F 2 = F3 = 2/3 
(that is, the bound for an optimal symmetric 1^2 doner) one has that F™ ax = 1/3 > 0. This means that, even if 
two optimal clones have been produced, there still remains some quantum information to produce a non-trivial third 
clone. A similar situation occurs for the case of cloning with discrete variables, as pointed out in Ref. 29] . 

Similar results occur for the generic m case. In fact, it can be immediately shown by inspection that substituting 
rij = (to — l)/m (that is, the bound for the noise introduced by an optimal symmetric 1 — > to doner) in Eq. (|37|l one 
obtains n" lm = (m — l)/m. Hence optimal symmetric cloning is recovered. Similarly, substituting rij — (m— 2)/(m— 1) 
(that is, the bound for the noise introduced by an optimal symmetric 1 — ► (to — 1) doner) in Eq. 137JI one obtains 
n min _ ^ m _ i)/( m _ 2), from which a fidelity F™ ax = 2 T ^~ 2 3 > follows. This confirms that the production of 
(to — 1) optimal clones still leave some quantum information at disposal to produce an additional non-trivial clone. 
An explanation for this effect may be individuated recalling that for large to the optimal doner coincides with an 
optimal measurement on the original state followed by to reconstruction 13]. As a consequence, one may expect that 
the production (reconstruction) of (to — 1) optimal clones leaves information (i.e., the measurement result) for the 
reconstruction of other ones. 

A question strictly related to the one faced above, and probably more significant from an information distribution 
viewpoint, is the following. Suppose that one wants to distribute the information encoded in the original state by 
fixing the ratio between the noise that affects all the to clones, and not by fixing the fidelities of (to — 1) clones. More 



specifically, suppose that one wants to give the minimum noise to, say, the first clone (rij > n\ for every j = 2, 



,to). 



Now fix the noise that affects the other clones by fixing their ratio qj with respect to the first one, that is rij = qj n\ . 
Then, which is the minimum noise n™ 111 allowed by our protocol for fixed qj? Solving Eq. I|37|l for m, one may find 
the following closed expression for n™ ln as a function of qj : 



to — 1 



-(m-lMl+Ejfc 



(41) 



The state \$f m ) that provides this optimal result is simply given by setting the Ni and Nj'a obtained by substituting 
back nf n = qj nf n in Eq. @B) and Eq. " 

As a final remark we point out that a general bound for the fidelities in a fully asymmetric 1 — * m cloning of 
coherent states has not yet been derived when m > 3. As a consequence, we cannot judge if the telecloning process 
introduced above is in general optimal or not for to > 3. Nevertheless, there are valuable indications for its optimality, 
i.e. the fact that it is optimal in the case of m = 2, and, as we have already pointed out, it is optimal for any m in the 
symmetric case. In addition, as already mentioned, our telecloning protocol allows to built a non-trivial additional 
clone when (m — 1) optimal ones have been produced. 



V. TELECLONING IN A NOISY ENVIRONMENT 



The protocol described in the previous section is referred to the case of ideal generation and propagation of the 
states \Sf? m ) as well as to double-homodyne detection with unit quantum efficiency. In order to take into account 
the possible losses and noise in the various steps, it is useful to reformulate the whole protocol in the phase space. 
Consider the characteristic function associated to the states |\P m ): %[<T m ](A) = cxp{— ^A T cr m A}. The covariance 
matrix cr m given in Eq. © can be written accordingly to the following bipartite structure 

(42) 

where A is a 2 x 2 matrix corresponding to mode do, while B and C are 2to x 2m and 2 x 2to matrices respectively. 
Consider now a generic Gaussian POVM, acting on modes oq and 6, defined by a covariance matrix M and a vector of 
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first moments X, i.e., \[M , X](A) = exp{- 
introduced above, corresponds to 



2 1V 



MA—iA T X}. The case of the ideal double-homodyne measurement 



M 



X 



(43) 



where <x in and X are the covariance matrix and the vector of first moments of the input state (mode b), whereas 
Z = {Ste[z\, 3m[z]} is the measurement result [we recall that P = Diag(l,— 1)]. Then, the state conditioned to the 
result Y is given by a Gaussian state with covariance matrix 



<t c = B-C t (A + M)- 1 C 



(44) 



and vector of displacements H — C T (A + M)~ 1 X . The protocol is now completed with the proper generalized local 
displacement introduced in the previous section, i.e., U z = <^) /l=1 DJ^(z). Averaging over all the possible outcomes we 
finally obtain the following expression for the covariance matrix of the Gaussian state at the output [30l | : 



B 



\A + M)PJJ 



"C-C 1 



(45) 



where J is given by the 2 x 2m matrix JJ = (1, . . . , 1). 

As already pointed out in Sec. lIIII if we consider a realistic scenario for the application of the telecloning protocol, we 
must take into account that the generation and the propagation of the states \^f m ) are affected by thermal background 
and losses. In particular, concerning propagation we can consider that modes a±, . . . , a m propagate in noisy channels 
characterized by the same losses T c . We may then define an effective propagation time r c = T c t equal for all the 
clones, while the effective propagation time To = Tot for mode oo is left different from t c . Consider in fact a scenario 
in which one has two distant location (see Fig. the sending station, where the double-homodyne measurement 
is performed, and the receiving station, where the clones are eventually retrieved. The distance between the two 
stations can be viewed as a total effective propagation time t t which can be written as t t = tq + t c . Then, the choice 
made above corresponds to the possibility of choosing at will, for a given t t , which modes (ai, . . . , a m or ag) will be 
affected by the unavoidable noise that separates the sending and the receiving station and to which extent. With a 
slight abuse of language, we may say that one can choose whether to put the source of the entangled state |* m ) near 
the sending station (t t — t c ), near the receiving one (t t = tq), or somewhere in between. A similar strategy has been 
pursued in |3J to optimize the CV teleportation protocol in a noisy environment. In the following, we will see how 
to determine both the optimal location and the optimal \^ m ) for a given amount of noise. For the sake of simplicity, 
the thermal photons /j, will be taken equal in all the noisy channels. As it is natural to expect, in the generation 
process all the modes will be also considered to be affected by the same amount of noise, characterized by v mean 
thermal photons. As a consequence, the matrix a m in Eq. l|4*2*|l should be substituted by (see, e.g., Ref. 0) its noisy 
counterpart: 



a/2 



(46) 



where we have used Eq. (|12|) . and defined 



'1©L e 



= (p+ ±)l 2ri 



(47) 



Performing the calculation explicitly, upon defining 7 C = e Tc , 70 = e T ° , k = \i + \ and £ = 1 + 2v, we obtain: 



A 

C T 



(48) 



where A = C [70M) + f (1 - 7o)l], C = C ^ToT^C and 

^7cATi + f(l-7c)l 7cBi,2 

7c»i,a 7cAT 2 + f(l-7c)l 



B = C 



\ 



lcB m -X,m IcAfm + f (1 - 7c) 1 J 



(49) 



A non-unit efficiency rj in the detection stage corresponds to have the covariance matrix of the double-homodyne 
detection given by M = P<Ti n P + |A1 (where A = Finally, considering an initial coherent state j^] and 

recalling Eq. H45[). we have M — i(l + A)l, whereas the covariance matrix of the m output modes now reads: 



B 



\A + M)PJ 



F T PC - C T I 



(50) 
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which in turn gives the following covariance matrix for the h-th clone: 

'M=(^-s)l. (51) 
In the Equation above, Fh represents the fidelity between the h-ih clone and the original coherent state: 
F h = {Det [<r h , n + ±l]y 1/2 

= |l + y +2k + C 70 (n + i - ^ +7c (jVj, + \ ~ - 2 V7o7 c WhCNo + l) } • 

A. Optimization of the symmetric protocol 

In order to clarify the implication of the formula (|52|l . let us focus our attention to the case of symmetric cloning 
(recall that in this case Ni = . . . = N m = N). Upon defining x = 4 — |, j T — e~ TT = 7o7 c and the following function 

/(iV,7o;£,7r) = ^(N-x)+ l0 (mN-x)~2^/ lT N(mN + l) , (53) 
7o 



(52) 



the fidelity reads as follows 



r = {{f(N,~/ ;x, 7T )+2 K +l + ^\ . (54) 



Our aim is now to optimize, for a fixed amount of noise, the shared state |^ TO ) and the location of its source between 
the sending and the receiving station. Namely, one has to find N and 70 which maximize the fidelity F for j T , k, 
C, A fixed. This, in turn, means to minimize f(N, 70; x, 7 T ) for fixed j T and x. The domain where to perform the 
minimization is the region N > and 7t < 70 < 1- We will see that the possibility of varying 70 will reveal crucial 
in order to adapt the ideal cloning protocol, presented in Sec. II VI to a noisy environment. 
Calculating the stationary points of f(N, 70; x, j T ) one finds: 



I at - , iTX 

s\ = <N = — , 70 



1 — x(m — 1) ' \ x + I 



I , T x jr (1 + mx) 

si = < N = —r— — -r , 70 = \ } . (55) 

m[l + xym — 1)J V mx \ 

The points s x and s 2 belong to the domain for | 7t < ^ ,x < and for | 7t < ,Vx| respectively. By 

evaluating the Hessian matrix associated to /(-/V, 70; x, J T ), it follows that both s± and S2 are not extremal points. As 
a consequence one has to look for the minimum of f(N,"fo;x,"f T ) along the boundary of the minimization domain. 
Three local minima are found in the three regions parametrized by 70 = j T , 70 = 1 and N — > 00, whereas the forth 
extremum is a maximum. In particular, in the first region the minimum is attained for 

7t > 1/m 



N = { m \ 1 (56) 
7t < 1/m 



m{m^ T — 1) 

Concerning the second and the third region, one finds that the minima are located at 



and at 



N = lT (57) 

m(7 T — to) 



7o = J-, (58) 

V TO 
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respectively. By evaluating the value of f(N, 70; %, 7t) in the minima, one eventually attain the global maximum 
F max of the fidelity. A summary of the results is given in Tab.Q] where we have reintroduced the effective propagation 
times r T , ro and defined the following quantities: 

F * = 2 JH (59) 

-2-4v + m{A + 2[2 + fi+u+(u- ^i)e- T T]} ' v ' 

J F lb ^2{A + 2(2 + ^ + i/)-2(l+/i + ^)e- rT } -1 , (60) 

F c = i2 + j + 2n-^^[l + ^ + iy + m(p-iy)}\ . (61) 



X 




opt 


N opt 




Vx 


< t t < lnm 


T T 


1 

m(me~ T r — 1) 


pa 


~ < x < 


r T > lnm 


| (t t + In m) 


TV -> 00 


F c 




mm < t t < In = — 

ma; 2 


| (t t + In m) 


TV —> 00 


F c 


TT >i n ii±4! 


r r 


e ~T T 

1 — m e~ T T 


F b 


x > ^=1 


r T > lnm 


t t 


e -T T 

1 — me~ T T 


F b 



TABLE I: Values of the optimized N° pt and r opt for fixed values of r T and x. The value reached by the fidelity ^ max for these 
optimal choices is given in the last column. 

An inspection of Tab. [I] shows very interesting features of the telecloning protocol in presence of noise, ft is 
immediate to recognize that the optimal value A opt is significantly different from the optimal value in the ideal case 
[Eq. J32J]- As a matter of fact N opt is divergent in some cases. Remarkably, in the optimization of N and ro the 
homodyne detection efficiency A plays no role, whereas the thermal noises /x and v introduce a dependence on x, i.e. 
only on their ratio. Furthermore, one may note that what we have called the best location of the source (that is TQ Pt ) 
is never given by the simple choices tq = or tq = t t /2. In order to clarify this point let us first consider the case 
t = 0, which can be physically implemented by homodyning mode a immediately after the generation of |^ m ), and 
then letting the other modes propagate to the receiving station where they are eventually displaced. An immediate 
calculation shows that in this case the fidelity (|54|) is maximized for N° pt = l/m(me TT — 1) and is given by 

2 m 

pmax/ T _ q\ _ " (Q2) 

V0 ; 2 e T r [1 + m((i - v) + 2 v]-m[A + 2(2 + + ' V ' 

Concerning the case ro = r r /2, whose physical implementation simply means to put the source of |\& m ) in the middle 
of the transmission line, one has that the fidelity is maximized for 7V opt = l/m(m — 1) and reads 

V 777 

pmax/ _ ley\ _ ^ " b (f,^) 

{ T/ ' m(4 + A + 4 /i) - 2e _TT / 2 [1 + 2^ + 2 m(/i - v)] V ' 

Notice that only in this case the optimization over N leads to the same |\& m ) as in the ideal case (see Eq. (|33fl ). A 
comparison of the last two instances with the optimal one, shows how significantly the choice of ro affects the value of 
the clones' fidelity. In Figs.03E] a nd[Hlwe compared the two fidelities given in Eqs. (|62[) and (j63|l with the one given 
in Tab. Q] (see captions for details). We clearly see that the optimized fidelity is much larger then the other two, thus 
providing a cloning beyond the classical limit for higher propagation times r T . As it is apparent from Fig. [5] we have 
F b < h Vr T . Indeed, it can be shown analitically that F° < ^ in any regime for which F max — F b . 

Besides what we pointed out above, the most striking feature of the proposed telecloning protocol is that it saturates 
the bound for optimal cloning even in the presence of losses, for propagation times r T < lnm, hence divergent as the 
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1 2 3 4 5 

FIG. 3: Comparison of the fidelities given in Eqs. 16211 (dotted line) and 1631 (dashed line) with the one given in Tab. Q (solid 
line). As an example, we have chosen the following parameters: /j, — 0.03, v = 0.01 and A = 0.02 (x = 1/51). The plots 
are referred to the case m = 2, 5, 10 and the vertical line corresponds to t t = In 2, r T = In 5 and r T = In 10 respectively. 
Accordingly to Tab. HJ the optimal fidelity is given by Eq. 1591 and Eq. 1611 at the left and at the right of the vertical lines, 
respectively. 
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FIG. 4: Comparison of the fidelities given in Eqs. 1621 (dotted line) and 1631 (dashed line) with the one given in Tab. Q (solid 
line). As an example of the case \jl < v, we have chosen the following parameters: /i = 0.05, v = 0.2 and A = 0.05 (x = —3/28). 
The plots are referred to the case m = 2, 3, 5 and the vertical line corresponds to r T = In 2, r T = In 3 and r T = In 5 respectively. 
Accordingly to Tab. HJ the optimal fidelity is given by Eq. 1591 and Eq. 16H at the left and at the right of the vertical lines, 
respectively. 



number of modes increases. More specifically, consider the first row in Tab. [I] and set fx = v = A = 0. Then, one has 
that for t t < lnm the maximum fidelity is given by i^ max = m/(2m — 1). That is, the optimal fidelity for a symmetric 
cloning can still be attained, carefully choosing TV and To- Such a result cannot be achieved letting the input state 
propagate directly to the receiving station and then cloning it locally [30| . Thus, in the context of our protocol the 
entangled resource significantly enhances the capacity of distributing quantum information. This is due to the fact 
that the transmission through a lossy channel of an unknown coherent state irreversibly degrades the information 
encoded in it, thus avoiding the local construction of optimal clones at the receiving station. On the other hand, 
multimode entanglement is robust against this type of noise and, even if decreased along the transmission line, it is 
still sufficient to provide optimal cloning. Actually, there is no need of an infinite amount of entanglement to perform 
an optimal telccloning process. 

As concern the case of higher transmission times, i.e. t t > lnm, the fidelity reads (again for ji = v = A = 0) 

F max = ^2 - ]J^Pj • (64) 

Eq. (|64|l shows that the fidelity is always greater then the classical bound F = h , which in turn means that the state 
used to support the protocol is entangled for any t t . This is reminiscent of the result already pointed out in Sec. IIIII 
where the full inseparability has been proved for any r T for m = 2 (notice that tq — r T /2 in case of Sec. 1111(1 . Here, 
we proved that the same conclusion is valid for any m. 

Another interesting feature in the case r T < lnm is that F max does not depend on t t if \i = v. Indeed, it turns 
out that for /i = and v ^ it is better to let the entangled resource propagate (up to tt = In m) instead of using 
it immediately after the generation. This effect may be naively understood by considering that the entangled state 
generated for v ^ is mixed and, as consequence, the propagation in a purely dissipative environment acts like a sort 
of purification process on it. As it is apparent from Eq. I|59|) this effect is present whenever /i < v (see also Fig.0J. 

Finally, a comment is needed concerning the scaling of the fidelity with respect to the number of modes m. We 
have already pointed out that for the case /j, = v = A = the fidelity remains optimal for times t t diverging with 
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FIG. 5: Comparison of the fidelities given in Eqs. I|62[l (dotted line) and l|63p (dashed line) with the one given in Tab. □ (solid 
line). As an example, we have chosen the following parameters: n = 0.6, v = 0.01 and A = 0.1 [x = 59/102). The plots are 
referred to the case m — 2, 3 and the vertical line corresponds to t t = In 2 and tt = In 3 respectively. Accordingly to Tab.Q] 
for m = 2 the optimal fidelity is given by Eq. 1591 and Eq. 1(> 1 1 at the left and at the right of the vertical line, respectively. 
For m = 3, it is instead given by Eq. 16011 at the right of the vertical line. 



the number of modes. However, when thermal noise is added (p, i/,A^0) the fidelity goes below the classical value 
F = | for times t t that become smaller as m increases, as it is apparent from Figs. El El and [5J Indeed, this is 
consistent with the fact that the optimal fidelity itself approaches the classical value F — \ as m increases. Hence, 
even a small amount of thermal noise is enough to cancel the benefits due to quantum entanglement. 



VI. CONCLUSIONS 



In this paper we have dealt with the properties and applications of a class of multimode states of radiation, the 
coherent states of group SU(m, 1), which represent a potential resource for multiparty quantum communication, as 
recent theoretical and experimental investigation have shown. In particular, the common structure of these multimode 
states allowed to consider a 1 — > m telecloning scheme in which a generic coherent state of SU(m, 1) plays the role 
of entangled resource. Exploiting the possible asymmetry of SU(m, 1) coherent states we have suggested the first 
example, in the framework of CV systems, of a fully asymmetric 1 — ► m cloning and have found the optimal relation, 
within our scheme, between the different fidelities of the clones. In particular, we have shown that when (m — 1) 
optimal clones are produced (accordingly to the general bound imposed by quantum mechanics), there still remains 
some quantum information at disposal. In fact, our protocol is able to use the remaining information to realize a 
non-trivial m-th clone. Our asymmetric scheme is aimed at the distribution of quantum information among many 
parties [3^ , and may find application for quantum cryptographic purposes [34| . 

In view of possible applications of our protocol in realistic situations, we have considered the effects of noise in the 
various stages of the protocol, i.e. the presence of thermal photons in the generation process, thermal noise and losses 
during propagation, and non-unit efficiency in the detection. We have derived the fidelities of the clones as a function 
of the noise parameters, which in turn allowed for adaptive modification of the protocol to face the detrimental effects 
of noise. In particular, we have shown that the optimal entangled resource in the presence of noise is significantly 
different from the one in the ideal case. Also the location of the source plays a prominent role. In fact, we have 
demonstrated that the optimal location is neither in the middle between the sender and the receiver, nor at the sender 
station. A striking feature of the optimized protocol is that, even in the presence of losses along the propagation line, 
the clones' fidelity remains maximal, a result which is not achievable by means of direct transmission followed by local 
cloning. This happens for propagation times that diverge as the number of modes increases. We then conclude that 
our optimized telecloning protocol is robust against noise. 
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